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Abstract 



We examine several different types of five dimensional stationary space- 
times with bulk scalar fields and parallel 3-branes. We study different methods 
for avoiding the appearance of spacetime singularities in the bulk for mod- 
els with and without cosmological expansion. For non-expanding models, we 
demonstrate that in general the Randall-Sundrum warp factor is recovered in 
the asymptotic bulk region, although elsewhere the warping may be steeper 
than exponential. We show that nonsingular expanding models can be con- 
structed as long as the gradient of the bulk scalar field vanishes at zeros of 
the warp factor, which are then analogous to the particle horizons found in 
expanding models with a pure AdS bulk. Since the branes in these models 
are stabilized by bulk scalar fields, we expect there to be no linearly unstable 
radion modes. As an application, we find a specific class of expanding, station- 
ary solutions with no singularities in the bulk in which the four dimensional 
cosmological constant and mass hierarchy are naturally very small. 

I. INTRODUCTION 

The Randall-Sundrum (RS) model for a warped 5-D geometry can account for the 
mass hierarchy, while providing a fine-tuning mechanism for canceling the 4-D cosmological 
constant without requiring a vanishing 5-D vacuum energy . Variants of the model have 
been constructed in which the cosmological constant is exponentially small and in which 
both the cosmological constant and the hierarchy problems may be solved simultaneously 
0]. However, the RS model and similar models that solve the hierarchy problem are un- 
tenable: the radion-mediated interaction of matter on the visible brane dominates over the 
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gravitational interaction 0,^. This has the consequence that matter on the visible brane 
(with positive brane tension) gives a negative contribution to the square of the four dimen- 
sional Hubble parameter for homogeneous cosmological models . Another manifestation 
of the problem is that the spacetimes are dynamically unstable: the 3-branes which they 
contain, whose positions must be arranged carefully to reproduce the mass hierarchy, move 
away from those special positions when slightly perturbed HJ^. 

As suggested by Goldberger and Wise (GW) 0, the introduction of bulk scalar fields 
can solve these problems by giving a mass to the radion mode and stabilizing the positions 
of the branes. This is consistent with phenomenology if the radion mass is at or above 



the electroweak scale []T0|. DeWolfe, Freedman, Gubser and Karch (DFGK) [|Tl|] showed 
how fully consistent stationary spacetimes could be constructed with bulk scalar fields, 
including the gravitational back-reaction of those fields; there have been further studies of 
spacetimes with bulk scalars along similar lines |lT2|jr3| . DFGK argued that the requirement 



that the induced metric on the brane be fiat requires one fine-tuning of the parameters 
describing the solutions. Generalizing the class of solutions by allowing a non-zero effective 



4-D cosmological constant A4 relaxes this fine-tuning |TT 

The purpose of this paper is to extend the work of DFGK [|ll|] and Kakushadze [|13| , 
to discuss some general properties of 5-D spacetimes with a bulk scalar field, and to apply 
those properties to the solution of the cosmological constant and the hierarchy problems. 
The main points of this paper are as follows: 

• The principal problem associated with the introduction of a bulk scalar field is that 
generic stationary solutions contain timelike curvature singularities in the bulk at fi- 
nite distances from the branes. This occurs both for static solutions (A4 = 0) and 
for stationary solutions with cosmological expansion (A4 7^ 0). Such singularities have 
already been encountered in work on self-tuning of the cosmological constant [1^. It 
is possible to avoid these singularities in two ways. One well-known way is to simply 
orbifold or otherwise compactify the fifth dimension in such a way that the singularity 
is never encountered. A second way, which is perhaps preferable, is to carefully choose 
the scalar field potential in such a way that the occurrence of singularities is prevented. 
In Sec. H below, we discuss a simple condition (first considered by Kakushadze [p!3| ) 



for identifying such preferred potentials for static spacetimes, and present several ex- 
amples. One interesting example is obtained by compactifying a 11-D spacetime down 
to 5-D with a single radion field which acts as a bulk scalar, as well as a 7-form field 
strength in 11-D descends to a non-dynamical 5-form field strength in 5-D that gener- 
ates a potential for the scalar field (Sec. |III C 1| below). By adjusting the 5-form field 
strength, this model can be rendered nonsingular. 

In the RS model, the metric's "warp factor" falls off exponentially as one moves away 
from the Planck brane, and this exponential fall-off underlies the RS solution of the 
hierarchy problem. In many models with scalar fields, the fall-off of the warp factor 
is faster than exponential, which facilitates solving the hierarchy problem (see the 
models in Sees. [Ill A| and [llllj| below). This point was mentioned in passing for a 
specific model by DFGK. 

We show that allowing the four dimensional cosmological constant A4 to be non-zero 
exacerbates the tendency to form spacetime singularities in the bulk. In particular, 
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models with bulk scalar fields that are nonsingular when A4 = should be expected 
to become singular for nonzero A4, although the corresponding singularities are rather 
mild. The singularities are weak enough (see Eq. |^0|) to be removed by simply re- 
quiring the bulk to become effectively AdS as the warp factor A{y) 0. In Sec. |^ 
below we derive a general method of constructing nonsingular models with cosmolog- 
ical expansion based on this idea. 

We can now construct models without bulk curvature singularities. With us living on 
the visible brane, it is possible to account for the extreme smallness of the electroweak 
scale and the cosmological constant simultaneously (Sec.0 below). These models are 
analogous to those of Refs. with the additional feature that the stability problems 
of Refs. BHI have been cured. 



II. SETUP AND GENERAL CONSIDERATIONS 

In this section we outline the general framework, review the procedure introduced by 
DFGK for constructing static and stationary solutions, and derive criterion under which 
bulk singularities do not occur. 

A. Basic equations 



We consider 5-D gravity plus a bulk scalar field with parallel 3-branes, for which the 
action is 



S 



J d'^xdyy 



1 p (V0)^ ... 



b •'yb 



(1) 



Here the coordinates are {x'^,y) for < /i < 3, the bth brane is located at y = yt, gab 
is the 5-D metric and g'^^^^ is the induced metric on the 6th brane. The brane tensions 
and potential V are functions of the bulk scalar 0, and is the 5-D gravitational coupling 
constant. 

We seek solutions to the field equations of the form 



ds^ = dy^ + A{y) -dr + e^"^6ijdx'dx^ 



(2) 
(3) 



although we shall concentrate initially on the static case H = 0. For this ansatz, the Einstein 
and scalar field equations reduce to |[Tl[ 



u' = 


2/«2((/,')2 


2H^ 


3 


A 


u^ = 




3 


3 


" + 2m0' = 


dV{ct^) 


dab 
_d(j)_ 



J2^b{(l))S{y - yb), 



+ 



A ' 

Ky-Vb), 



(4) 
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where 

= ^ (5) 

and primes denote derivatives with respect to y. Only two of these three equations are 
independent: the scalar wave equation follows from the other two via the Bianchi identities. 
The jump conditions at the branes are 

u\f = Wm u - Wm u = -2g,, = a'M, (6) 

y-*yb+ y-^Vb- Vb 

where 0b = (t){yb) is the value of the bulk scalar on the 6th brane and 

Qb = ^KVb(0b). (7) 



The Ricci scalar is 



R = -Au' - + 1^ (8) 
^.W + ^^^. (9) 



B. Method of obtaining solutions 

We now review the method of generating solutions introduced by DFGK. For any given 
solution of Eqs. @, we can imagine inverting the relation = to obtain y as a function 
of 0. Leaving aside questions of single-valuedness, we can imagine changing independent 
variables from y to 0. Since 0'(?/) is also a function of y, we can likewise take 0' to be a 
function of 0. By analogy with the practice in supergravity theory, let us define a function 
by 

^' - (10) 

^ 2 d(h ^ ' 



In the H = case, the first of Eqs. (§) then implies that, away from branes, 

90 ~ 3"*^ ~ ~y d(j) ' 

so that 



(11) 



u = -K^W{(l))/3. (12) 
The second of Eqs. (|) then becomes 



l^'w{<pr. (13) 



The procedure for obtaining solutions in the H = case is (i) choose a potential V{(j)) and 
superpotential W{(j)) that are related by Eq. (|T3p; (ii) integrate Eq. ( pUD to obtain as 
a function of y; (iii) combine this with Eq. (p!2|) to obtain m as a function of y; and (iv) 
integrate Eq. to obtain the warp factor A as a function of y. 
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C. Occurrence of singularities in the bulk 



Combining the definition (|TUp of tfie superpotential W witfi tlie expression (|^) for tlie 
Ricci scalar gives 



(14) 



Eq. ([T^) is applicable to models with H 0. When H = 0, substituting Eq. ([T3| ) into Eq. 
(0) yields 



2k' 



2 



5/s:W^(0) 
9 



(15) 



Eqs. ([I^ ) and (|15D show why singularities are rampant in models with bulk scalar fields. 
For many choices of W{(j)) and V{(j)), it turns out that |-R(0)| — > oo as |0| oo. Examples 
include V{^) = /iW2 or = Xcp'^/A or 1/(0) = Voe''^/^^, all of which are simple, 

and even well-motivated in some respects. Thus, unless (p is constrained to remain finite 
everywhere, a singularity will be encountered. 

To see how to avoid singularities in general [|T^, consider first a spacetime with a single 
brane, a.t y = i/b, and suppose that the value of the scalar field on that brane is (pb- From 
Eq. (p!0D, we can solve for ?/ as a function of 0: 

y(0) =y, + J(0b,0), (16) 

where 

Now if there exist finite values 0_oo and 0+oo of the bulk scalar field such that \I{4>b, (f>±oo)\ = 
oo, then the spacetime can be infinite in the fifth dimension with — >■ 0±oo as y —oo 
and ?/ — >■ oo. In such cases, an infinite range in y is mapped onto a finite range in 0, thus 
potentially avoiding any singularities. 

From Eq. (|T7|) , divergences of 1(06,0) will occur only at zeros of W {(/)). Suppose now 
that W'{(j)) has a zeros at both > 0b and < 0b, and that W'{(f)) vanishes at those zeros 
at least linearly (but not like |0 — 0±oo|^ with p < 1). If, in addition, neither W'{(f)) nor 
\/(0) diverges in between these zeros, then the spacetime can extend infinitely in the fifth 
dimension away from the brane in either direction, without any singularities. A similar 
result can be obtained for spacetimes with more than one brane. In this case, the condition 
is that I{(j)bi<Pb') be finite for any two branes 6,6', and that for the two bounding branes, 
6±, there exist 0±oo such that I{(f)b±, 0±oo) = ±oo with no intervening infinities in V{(j)) and 

91^(0)/90. 

In models where these conditions for avoiding singularities are met, V{(j)) tends to a 
constant asymptotically, as a zero is approached. Thus, such models, for H = 0, tend toward 
the RS solution far from branes, provided that the scalar field potential is negative. However, 
there will also be regions in such solutions where the spacetime metric differs considerably 
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from the RS model, and where the warp factor drops far faster than exponentially. We 
shall show how this comes about in a specific model below. The H ^ solutions may 
also tend asymptotically toward their counterparts with uniform bulk cosmological constant 
However, the relatively harmless particle horizons characteristic of those uniform bulk 
cosmological constant models will generically be transformed in to curvature singularities, 
unless there are fortuitous cancellations, as can be seen from the A limit of Eq. (|^). 

Models where singularities would be inevitable in an uncompactified geometry can be 
truncated by compactification to a region that does not include any singularities. Such 
models need not approach the RS model, even asymptotically, and can have warp factors that 
vary much more rapidly than exponentially throughout the compactified fifth dimension. 



III. MODELS WITHOUT COSMOLOGICAL EXPANSION 

In this section we discuss several different choices of superpotential W{(j)) and the prop- 
erties of the corresponding solutions in the static case H = 0. 



A. Even superpotential 

The RS model corresponds to W{(j)) = ±a, where a > is a constant; this model has 
V = — DFGK showed that the GW stabilization mechanism can be implemented if 
instead W{(f)) = ±(a — 60^), where 6 > is another constant, for which 

vi^^^-^jf^*i(,^-'-^y^. (18) 

Choosing the + sign in W{(f)) for y > (and vice- versa), DFGK found (f){y) = 0(O)e~''^, and 
therefore (via Eq. [Q) Su/k^ = —a + 60^(O)e~^^^, which implies a transition between two 
regions of constant u, one at small by and the other at large by, where the model becomes 
equivalent to the RS model. This model has dW {(p) / d(j) = at = 0, and is nonsingular 
as long as the maximum value of is finite. Thus, it avoids singularities by mapping an 
infinite range of y to a finite range in (p, as was discussed in Sec. [11 G. 



DFGK noted that they could also choose W{(f)) = ±(a + 60 ), with 6 > 0, for which 



VW-~-^-^{^'-^-^]--^. (19) 



Eqs. ([TSp and (|T9D can be rewritten in the form 

VW = -A,^t^-^, (20) 
where the =F signs apply to Eqs.(|l8D and ([T9| ) respectively, and where 



6 



Despite appearances, the differences between the two models ([T8|) and (^) are not minor. 
For the + model, we have W{(f)) = a + hcl)^ for ?/ > 0, so that 0' = 60, 0(?/) = 0(O)e*^ 
and ?)u/ = —a — 60^(O)e^''^. This yields a warp factor A{y) which is an exponential of an 
exponential: 



1 o 1 



In A{y) = In A{0) - -n'ay - -fi:^0(O)' [e'"^ - ij . (22) 

This super-warped version of the GW model, obtained from an apparently insignificant 
change to the potential (^0]) , becomes singular if the fifth dimension is not compactified, 
because R{4>) ^ — oo as ^ oo. Any compactification of the fifth dimension avoids the 
singularity by construction. The model can be altered slightly to become bulletproof against 
singularities by taking W{(j)) = ±(a + 60^ — 60'*/20^). For this choice of W we have 

so that 0OO and u — (K^/3)(a + 60^/2) at large y, reducing asymptotically to the RS 
model. Note that the total change in the logarithm of the warp factor in the super-warped 
regime is limited to ~ ^^0^; which may be large. 



B. Odd superpotential and Gaussian warp factor 



The DFGK realizations of the GW stabilization mechanism for the RS model are based 
on choices of superpotential W{(f)) that are even functions of 0. There are also simple models 
in which W{(j)) is odd, first considered in Ref the simplest is 



W{4>) = ±260, 

with 6 > 0. (A constant term in W{(j)) can be absorbed into the definition of 
model (121), we find that 



(24) 
For the 



u 



In A{y) 



m + by, 

2^260(0) 
3 

InA(O) 



2K'b'y 



2K%(j){0)y K%'^y'^ 



(25) 



for y > Q. Thus, the warp factor in this model is Gaussian rather than exponential, and 
therefore falls off more rapidly than in the RS model. Note that A{y) has a maximum value, 
and all positive tension branes must be located away from that maximum. 

This model contains a singularity at y ^ oo since (f){y) diverges there. The singularity can 
be avoided by compactification, or can be removed altogether by setting W{(f)) = ±26(0 — 



/30^), in which case 



4>{y) = 0OO tanh 



tanh 



'm\ , by 



(26) 
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which tends to 0oo asymptotically. The warp factor for this model is given by 



lnA{y) = InA(O) - 



4^ 
9 



In 



cosh 



by 



1 



by 



[cosh.{by / (poo) + 0(0) sinh(6?//0oo)/0c 



(27) 



and tends toward In A{y) ^ — 4k^60oo1//9 for by /(poo ^ 1, which is the same scaling as in 
the RS model. The maximum change in the logarithm of the warp factor in the Gaussian- 
warped regime is of order 
above. 



similar to what we found for the model of Sec. pi A 



C. Exponential potential 

The models considered so far become singular only in the asymptotic regime y oo. 
Thus, whether or not these singularities merit a potential-altering cure is largely a matter 
of taste, since compactifying to any finite size in the fifth dimension removes the singularity 
without any sort of fine-tuning. However, there are other models which are "spontaneously 
singular" in the sense that they develop singularities at finite y. An example is when the 
superpotential is an exponential. 



W(0) = 2ae" 



(28) 



For this model, we have 



-k<f> 



gfc</.(0) _ ^2 



ay 



(29) 



which diverges aty = [k'^a) ^e^'^'^'^\ implying a divergence of the Ricci scalar. The singularity 
is removed by the following slight modification of the superpotential, 



W{(p) = 2a 



k 
— ( 
Q 



iq-k)4>oo-q4> 



where g > A; is a constant. The associated potential is 



V{<t>) 



-2fc((/--0oo) [ ^2 _ 



j^^-2q{4>-4>^) / ^2 _ 



21,2^ 



An^k 



For this model the relation between (p and y is given by 



k{4>-<Poa) 



rf 



.q/k-l _ I 



k^ae-'"^- 



(30) 



(31) 



(32) 



which can be integrated numerically in general. The integral can be done analytically in 
special cases. For example, for q = 2k, q = 3k, and q = 5k we find the results 
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gfc(<^-</>oo) _j_ ^k[(t>-4>oo) _ ]^ _ gA:(9i(0)-</)oo) _|_ Jj^ gA:(</)(0)-9ioc ) _ ^ — k"^ ae~'^'^°° y , 



and 



1 k, 
tann — 
2^ 



4 



gfc(0(o)-0oo) _^ _ In 



tanh-(0(O) 



tann — uz> — 



1 -1 / 

tan f 

9 V 



fc{<ii-</'oo) 



gfc(</.-</>oo) + _intanh 
4 



1 -1 / 
- tan I 



2 

M't'-'Poc) 



respectively. For large values of y this model reduces to the RS model, since 
above and Su/k,"^ — > — 2ae~'^'^°°(l — k/q), a constant, 



(33) 
from 



1. Specific realization obtained from dimensional reduction 



A potential reminiscent of the potential (^TJ) arises from compactification of an 11-D 
spacetime to 5-D; the result is 

V{^) = -A,e-'^/' - cge-^ + ye'^^. (34) 

Here it is assumed that the six compactified dimensions have a single associated radion field, 
ip, the curvature associated with the compactified dimensions is Cg > 0, A;, is the descendant 
of an 11-D cosmological constant, and £" is a 5-form field strength (descended from a 7- form 
field strength in 11-D). Although the potentials (^) and ( |3^ are similar, there are no choices 
of parameter values for which the potentials match up term by term. This is because the 
exponents in three terms in Eq. (R3) are in the ratios 1 : 1.5 : 3, whereas the exponents in 



Eq. (|3TD are in the ratio 1 : (A; + q)/{2k) : q/k. However, we can find parameter choices for 
which the potentials match if we allow one of the two parameters Cg and A;, to vanish. 

If we drop the term proportional to Af, and retain curvature in the six compactified 
dimensions, we find two possible identifications: 

q = 2k, fc2 = 2/cV3, ^' = C6 = (K'aV3)e-2'=<^-, (35) 

and 

q = 3k, = 4^727, = cg = (16/t^aV27)e-2'^°°. (36) 

If we set Cq = but retain A;, 7^ 0, there are also two possibilities: 

q = 3k, e=AK^/Q, A,= (4/€V/9)e-2^^°°, 8^ = {^^0" /27)e-^^^^ , (37) 

and 

q = hk, e=AK^/7h, Ab = (12/€2aV25)e^2fe<^-, S"" = {?>2^a^ /7h)e~^^'^^ . (38) 

Whichever possibility we choose, we can regard the combination t^?a'^e~'^^'^°^ as a derived 
quantity, determined by either A;, or ce, depending on which is nonzero. The 5-form field 
strength which involves a constant of integration, can then be adjusted to produce a non- 
singular model. The 5-form fields therefore may play a central role in eliminating spacetime 
singularities from 5-D models derived by dimensional reduction from 11-D. 
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IV. MODELS WITH COSMOLOGICAL EXPANSION 



A. Bulk singularities at zeros of the warp factor 

When the bulk geometry is anti-deSitter (AdS) space, metrics of the form (0) have non- 
singular particle horizons at zeros of the warp factor A{y). From the expression for the 
Ricci scalar, one might anticipate a singularity from the term ocl/A, but that divergence is 
cancelled by divergences in the u"^ and u' terms for an AdS bulk. 

However, when scalar fields are included, the surfaces where A{y) =0 generically corre- 
spond to curvature singularities. This is the case even when we set (j)' = ^dW{(j))/d(f) and 
choose dW{(j))/d(f) so that (f) remains finite over the entire range of y. The simplest way to 
see this is to consider the first of Eqs. @ in a neighborhood of a point where A{y) = 0. 
Suppose that (p' is finite and nonzero in that neighborhood, so that as long as we restrict 
attention to a small enough region, we can take it to be a constant. If = k^(0')^/3, then 

This equation would be exact for dW{(f))/d(j) = constant, which led us to the Gaussian model 
in Sec. |T^. Multiply by m = A' /A and integrate to find 

«2 =(^^J = -V In I A| + ^ + 4A;2. (40) 

Here fc^ is a constant which may be positive or negative, although we shall take it to be 
positive without exception to recover the RS model when q = = H. Substituting these 
results for u' and u"^ into Eq. (P), we find 

R = -20P + + 20g2 In \A\ (41) 

as y4 ^ 0. Thus, the Ricci scalar diverges at this point, which is a true singularity, unless 
= 0. The divergence is relatively mild (logarithmic). In fact, for H = this is the same 
divergence as was found asymptotically in the Gaussian model of Sec. pTl] . 

B. Method of constructing singularity-free solutions 

In the singularity-free models with if = 0, we saw that obtaining a mapping of the infinite 
y domain to a finite range of required that 0' — > asymptotically. For H = models the 
singularity is absent, because A ^ only asymptotically, and q^lnA — > asymptotically 
as well. For H ^ 0, nonsingular models can be constructed parametrically by the following 
procedure. Define a function g{A) by 

2k2(0')2 dg(A) , , 

= 2A^ > 0. (42) 

with (7(0) = 0. Then Eq. ( ppj ) generalizes to 
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u 




-AgiA) + ^ + 4fc2 



(43) 



Comparing this with the second of Eqs. (^ we obtain 

3 " dA + ~ 



(44) 



The procedure for obtaining a solution can now be summarized as follows: (i) pick a function 
g{A) and the integration constant fc^; (ii) solve Eq. (^) to obtain A as a function of y\ (iii) 
solve Eq. (|4^ ) to obtain as a function of (iv) use Eq. (^if ) to compute the potential 
1^(0) corresponding to the solution just obtained. By construction, the Ricci scalar 



R = -20k^ + 8A^ + 20g{A) 

is finite at A 0. 

Let us illustrate this procedure with the particular example 

g{A) = q'A^, 



(45) 



(46) 



with p > 0. (For small p, we expect this model to resemble the singular Gaussian model.) 
For the choice (^) we find 



(47) 



Apparently, there is a maximum value of A{y) for such models; for small H'^/k'^ we find 
Amax — {k/qY^'P + H'^/pk'^. Note that in this particular model, freedom to rescale A{y) 
also implies that q may be rescaled arbitrarily; this need not be true for g{A) that are not 
scale-free. When H = 0, the solution to Eq. (47) is (apart from an overall ambiguity in the 
sign of 0' and hence 0) 



k 



2/p 



gcosh pky I 



^ = tan- (e^'y) - ^ 



2k'^V{<P) 



Ak' 



-l+(l + |]cos2 




(48) 
(49) 
(50) 

= at y = 0. As y 

^k'^p/3 ranges from to Tr/2. Note that the maximum value of 
V{(f)) is positive in this model, and occurs at the maximum of A{y); this is already apparent 
from Eqs. (|43|) and (|4^) . In a particular realization of this model, with branes at specific 
values of 0, the maximum value of the warp factor may never be encountered. Finally, notice 
that for small values of p, the warp factor A{y) ~ e~^'^ ^ for pky ^ 1, and, for any value 
of p, A{y) ~ e~''y for pky ^ 1. These features are reminiscent of the modified nonsingular 



where we have located the maximum of A{y) at y 
ranges from zero to infinity. 



and chosen 
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Gaussian model of Sec. |T|, and the potential (|50|) therefore represents an alternative way 
of regularizing the Gaussian model to avoid singularities. 

When H = 0, Eq. (|D shows that A{y) only when y ^ oo. When H ^ 0, Eq. (g^) 
shows that l^'/^l is larger than for H = 0, and we therefore expect the position of the zero 
of A{y) to move inward from ?/ = oo to finite y. If yo is the zero, so that A{yo) = 0, then 
A[y) may be found by inverting the equation 

w N /"^ dA 

mvo -y)= I 51 

and the scalar field may be found from 

\AkV Jo k^-AP + H^A 



Here 0(0) is the value at A = 0, and we have assumed 0' > 0. Eqs. (|52|) and (|^) may 
be used to evaluate the potential V{(j)), which will have a nontrivial dependence on H. For 
A < H^/(^+p\ the scalar field plays an insignificant role, and we recover the vacuum solution 
A{y) ~ {H^/k'^) sinh2[fc(|/ - yo)]. For A > ij2/(i+p)^ ^j^g Q^j^g^ hand, we can neglect H , 
and we recover the solution (^) for A{y). Corrections to these approximate solutions may 
be obtained from Eq. (|5ll) . In general, the behavior of models with scalar fields is subtle for 
small values of H^, because nonsingular models require that the effects of the scalar field 
become insignificant as A{y) ^ 0, so that the terms become important in these regions 
despite the smallness of H^. 

These nonsingular models appear to share common features with models developed pre- 
viously without bulk scalar fields particularly the occurrence of non-singular surfaces 
where A{y) = 0. Presumably, these surfaces are particle horizons similar to those which 
occur in the pure AdS case. 



V. COSMOLOGICAL CONSTANT AND HIERARCHY PROBLEMS 

To examine the implications of these models for the value of the cosmological constant. 



focus again on the model (|6) and rewrite Eq. (53) as 

First, let us consider the S'^/Z2 orbifolded model, with a brane at each end, where the warp 
factors at these branes are Ai, i = 1,2, and the brane tensions are CTj. Note that it is possible 
for both brane tensions to be positive, if one of them is located a.t y > yo and the other at 
y < yo- We define the quantities Qi and Q2 by 

Using the jump conditions (|^) along with Eq. (^3]), where the constant term in V{(j)) is 
adjustable (for example, via unimodular gravity), we find that 
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^2 ^ A,Ql - A^Ql 



Ai-A 



2 



A,A,{Ql-Q\) 
^ = M-A, ■ 

In the regime A2 ^ Ai, these formulae reduce to 

H'^A2{Ql-Ql), (56) 

where we have assumed that Ql ~ Q^. Note that Eq. (^) requires that (^2 > Ql- The 
result (0) implies that the four dimensional cosmological constant A4 oc H^/Ai is sup- 
pressed by the ratio A2/A1 of warp factors, which will be vanishingly small for suitable 
brane separations. 

Now suppose that our Universe lives on a test brane located somewhere between the 
branes 1 and 2. Denote by Au the warp factor on our brane, and let brane 1 be the Planck 
brane. Then, we find that the electroweak scale on our brane is given by m|;,y = m1,Au/Ai, 
where mp is the Planck scale. The expansion rate on our brane is = H'^/Au, and 
therefore 

Hi, _ AjA^iQl - Qj) A.A^iQl - Q|) _ ^^^^ 



"^Ew mlAfjlAi - A2) im?pAl 

Assuming that ~ mp, we see that Hlj/m\-^ can be exponentially small provided that 
Ajj ^ A1A2. Roughly speaking, this will be the case if the test brane on which our 
Universe resides is closer to the Planck brane, brane 1, than it is to brane 2. Turning on 
a small positive brane tension for the visible brane does not change the above qualitative 
features. Even for models in which singularities are not absolutely prevented, they can still 
be avoided in the orbifold case, which requires a negative tension brane at one fixed point. If 
one only wants to solve the mass hierarchy problem, this negative tension brane is the visible 
brane f^. If one also wants to solve the cosmological constant problem, the visible brane 
must be identified with a third brane between the orbifold boundaries, as in the nonsingular 
models discussed above. 

This result can be generalized to other H spacetimes with no bulk curvature sin- 
gularities, and remains valid for uncompactified, multi-brane models as well. It can also 
be generalized to cases where the branes are charged under some 4-form potentials [Q. In 
this case, the constant term (and/or other coefficients) in V{(j)) depends on the 5-form field 
strengths, and the background values of these field strengths are determined by satisfying 
the jump conditions at the branes. With the warp factor normalized to unity at the Planck 
brane, A4 is proportional to Ap, the warp factor at the nearest particle horizon. This factor 
can easily be exponentially small. Since the visible brane must be closer to the Planck brane 
than the particle horizon, the warp factor A^ at the visible brane can easily be exponentially 
small as well, providing a solution to the hierarchy problem 0, but A^ » Ap, as observed 
in nature. 
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VI. DISCUSSION 



In this paper, we have constructed several different models of 5-D brane worlds with 
bulk scalar fields. We gave a general procedure for constructing static {H = 0) models with 
no curvature singularities in the bulk. A specific model was constructed corresponding to 
compactifying from an 11-D theory with a 7-form field strength down to 5-D, where the 
7-form field descends to a 5-form field strength. By adjusting the 5-form field strength, we 
were able to render the model nonsingular. In general, we found that nonsingular models are 
more warped than the RS model on relatively small distance scales, but tend towards simple 
exponential warping on large scales. Thus, the RS behavior is robust in the asymptotic 
regime, although stabilization of branes via the introduction of bulk scalar fields may require 
the branes to reside where the fields produce substantial extra warping of spacetime. 

For expanding [H ^ 0) models, we found that singularities occur even for those scalar 
field potentials that do not give rise to singularities in the H = case. We therefore 
were compelled to explore a different method for finding nonsingular models with H 0. 
A subclass of these models incorporates a non-singular surface on which the warp factor 
A{y) vanishes, analogous to the particle horizons that occur in the pure AdS case ||^. One 
specific class of such models, which corresponds to a sinusoidally varying potential in the 
nonexpanding case, was treated in some detail. For this case we constructed an orbifolded 
model, with one bounding brane being the Planck brane, and the other bounding brane (on 
which the warp factor is exponentially smaller) having either positive or negative tension. 
In this model, we live on a test brane between the two boundaries branes, and the expansion 
rate on our brane is exponentially smaller than the local electroweak scale provided that our 
brane is located somewhat closer to the Planck brane than to the other brane. This feature 
does not appear to depend in an essential way on details of our particular model, and even 
can be shown to hold when singularities are avoided (by orbifolding) rather than prevented 
outright by choice of potential. Thus, there are models in which the cosmological constant 
and mass hierarchy problems can be explained simultaneously. 

Two fine-tunings of the parameters of the RS orbifold model are required to get flat 4-D 
spacetime [Q. Allowing non-zero A4 relaxes one of these fine-tunings, while the other can be 
eliminated by allowing the bulk cosmological constant to depend on 5-form field strengths, 
whose piecewise constant values are determined by boundary conditions at the branes. When 
the branes are relatively far apart, A4 turns out to be exponentially small However, the 
models of Ref. @] were dynamically unstable. Bulk scalar fields can be introduced to fix 
brane positions. When A4 = 0, there are two sources of fine tuning in models with bulk 
scalar fields, one associated with expressing the potential V{(j)) in terms of a superpotential, 
W{(j)) [Eqs. ( pUj ) and (P^], and the other associated with adjusting W{(j)) so that the solution 
becomes nonsingular. Given a potential V{(f)), there is no guarantee that a W{(f)) can be 
found that satisfies Eq. (|T3p without tuning the parameters [^. For example, a general 



quartic potential of the form V{(f)) = ao+a20^/2+a40^/4 contains three parameters, whereas 
Eq. (|T^) or (^) contains only two, implying that a relationship among Oq, 02 and 04 is needed 
in order for V{(j)) to arise from a superpotential W{(f)). Moreover, Eq. (^) corresponds to a 
nonsingular model, whereas Eq. ([19|) is singular; the two independent parameters in a quartic 
potential have to be specifically adjusted to prevent singularities. The adjustment can be 
achieved by the introduction of adjustable 5-form field strengths, and we may speculate 
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that, in an evolving bulk, the 5-form field strengths may relax naturally (e.g. via bubble 
nucleation) to values that prevent the occurrence of singularities. (A specific example of 
the role that 5-form field strengths might play is given in Eqs. (|35D-(|38D of Sec. [Ill C 1| .) In 
general, additional fine-tunings required to prevent singularities need not be invoked in an 
orbifold model, where compactification may simply avoid the presence of singularities that 
could otherwise arise in an uncompactified model. 

When A4 7^ 0, some fine tuning is still needed, apparently, for singularities to be absent, 
but we had to resort to a different method, not based on a superpotential, to construct 
nonsingular solutions in this case. The problem arises because singularities can only be 
prevented if (p' when A{y) ^ 0, a condition that does not follow readily from a 
description in terms of a superpotential, and therefore represents a different kind of tuning. 
This suggests that the bulk alone may be supersymmetric, making a description in terms of 
a superpotential W{(f)) natural, but that the presence of branes breaks the supersymmetry 
completely, and leads to nonzero A4. As in the nonexpanding case, it is also possible to 
preserve models based on a superpotential without any special tuning of the parameters in 
W{(j)) in a appropriately compactified or orbifolded model, but it is never possible to relate 
V{(f)) directly to W^cj)) when A4 ^ 0. 
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